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Abstract 

An analysis of the non-modal growth of 
perturbations in a boundary layer in the presence of a 
streamwise pressure gradient is presented. The 
analysis is based on PSE equations for an 
incompressible fluid. Examples with Falkner-Skan 
profiles indicate that a favorable pressure gradient 
decreases the non-modal growth while an 
unfavorable pressure gradient leads to an increase of 
the amplification. It is suggested that the transient 
growth mechanism be utilized to choose optimal 
parameters of tripping elements on a low-pressure 
turbine (LPT) airfoil. As an example, a boundary- 
layer flow with a streamwise pressure gradient 
corresponding to the pressure distribution over a LPT 
airfoil is considered. It is shown that there is an 
optimal spacing of the tripping elements and that the 
transient growth effect depends on the starting point. 
The amplification is found to be small at the LPT’s 
very low Reynolds numbers, but there is a possibility 
to enhance the transient energy growth by means of 
wall cooling. 

Nomenclature 

E = energy norm 

G = E out /E in energy ratio 


H 

= shape factor 

h l 

= f’L/U eL 

h Ls 

= f L s /U exit 

H ref 

= scale in y- and z- directions 

L 

= length along the streamwise 
direction 

L ref 

= reference length 


Re 
R H 


T„, / T 


ad 


u 


eL 


ref 


U 

V 

W 

X 

y 

z 

p 

Ph 

s 


V 

e 

p 


Superscripts 


surface length 
Reynolds number 
U eL L/v 

temperature factor (ratio of the wall 

temperature to the temperature of 
the adiabatic wall) 
streamwise velocity component of 
the mean-flow velocity 

free-stream velocity at x = L 

reference velocity 

normal-to-the wall velocity 
component of the mean-flow 
velocity 

streamwise velocity disturbance 
normal velocity disturbance 
spanwise velocity disturbance 
streamwise coordinate 
coordinate normal to the wall 
spanwise coordinate 
spanwise wave number 

Hartree parameter 

f'aa ! E ref L ref 

kinematic viscosity 
momentum thickness 
pressure disturbance 


T = transposed 

Subscripts 

e = free stream 
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exit 

= exit conditions 

in 

= starting point 

out 

= ending point 


Introduction 

Laminar-turbulent transition in shear flows is 
still an enigma in the area of fluid mechanics. The 
conventional explanation of the phenomenon is based 
on the instability of the shear flow with respect to 
infinitesimal disturbances. The conventional 
hydrodynamic stability theory deals with the analysis 
of normal modes that might be unstable. The latter 
circumstance is accompanied by an exponential 
growth of the disturbances that might lead to laminar- 
turbulent transition. Nevertheless, in many cases, the 
transition scenario bypasses the stage of the 
exponential growth associated with the normal 
modes. This type of transition is called bypass 
transition. Observations of laminar-turbulent 
transition in plane Couette flow, in circular pipe flow, 
in boundary layers at relatively low Reynolds 
numbers, etc., can serve as examples of bypass 
transition. An understanding of the phenomenon has 
eluded us to this day. One possibility is that bypass 
transition is associated with so-called algebraic 
growth of disturbances in shear flows 1-3 . 

The phenomenon of algebraic growth of 
disturbances in shear flows has been of great interest 
during the last two decades as it may be associated 
with the bypass transition mechanism. Ellingsen and 
Palm 1 considered, in the inviscid case, an initial 
disturbance independent of the streamwise coordinate 
and found that the streamwise disturbance amplitude 
may grow in time, even though the basic flow does 
not posses an inflection point. Landahl 2 showed that 
all parallel inviscid shear flows are unstable to a wide 
class of three-dimensional disturbances. The result is 
independent of whether or not the shear flow is 
unstable to an exponential growth of wavelike 
disturbances. This type of instability that is not 
related to exponential growth is also referred to as 
“ non-modal growth. ” Mathematically, the effect of 
non-modal growth is associated with non-normality 
of the linearized Navier-Stokes operator and non- 
orthogonality of the eigenfunctions. In simple words, 
the essence of non-modal growth is the possibility of 
combining the exponentially decaying modes in such 
way that their sum will possess transient growth. One 
can find a vast bibliography on the topic in the Otto 
Laporte Award Lecture by Reshotko 3 and in a 
monograph by Schmid and Henningson. 4 

Numerical analysis of spatial non-modal growth 
within the scope of the linearized boundary-layer 


equations for an incompressible flow over a flat plate 
was carried out in Refs. 5 and 6. Spatial analysis 
within the scope of the linearized Navier-Stokes 
equations (quasi-parallel approximation of 
compressible and incompressible flows) was 
presented in Refs. 7 to 9. The main results of these 
theoretical models are as follows: 

• A system of counter-rotating streamwise 
vortices, which are periodic in the 
spanwise direction, provides the strongest 
growth of the disturbance. 

• There is an optimal spacing of the 
streamwise vortices, leading to the 
strongest effect. 

The effect of pressure gradients on the transient 
growth mechanism was considered within the scope 
of temporal theory by Corbett and Bottaro 10 and 
within the scope of spatial theory by Tumin and 
Reshotko. 9 These results were based on the quasi- 
parallel flow assumption. Tumin 11 analyzed the 
pressure gradient effect for the Falkner-Skan profile 
within the scope of an analytical model when the 
spanwise wave number is very small. The pressure- 
gradient effect within the scope of spatial theory with 
nonparallel base flow and finite spanwise wave 
numbers has not been considered, yet. 

Another motivation for the present work stems 
from separation flow control on low-pressure 
turbines (LPT). The performance of LPTs is strongly 
affected by the flow separation. There is a possibility 
of delaying the boundary-layer separation by tripping 
the boundary layer with the help of roughness 
elements or other devices. Usually, trial-and-error 
method is used to determine an appropriate 
placement of the control elements. This approach is 
time consuming and expensive. A recent 
investigation by Reshotko and Tumin 12 demonstrated 
that roughness-induced transition might be related to 
the transient growth mechanism. 

Periodically spaced in the spanwise direction, 
roughness elements generate a system of counter- 
rotating streamwise vortices. Due to a secondary 
instability mechanism, the streamwise vortices can 
lead to earlier transition to turbulence. They also 
provide a mixing enhancement due to redistribution 
of the streamwise momentum. Consequently, 
optimization of the streamwise vortices for maximum 
energy growth leads to maximizing of the flow 
control effectiveness. In the present work, analysis of 
the optimal disturbances/ streamwise vortices 
associated with the transient growth mechanism will 
be performed for boundary layers in the presence of a 
streamwise pressure gradient. The theory will provide 
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the optimal spacing of the control elements in the 
spanwise direction and their placement in the 
streamwise direction. 


Governing Equations 

Because the flows of interest have relatively low 
Mach numbers, we consider steady three-dimensional 
disturbances in an incompressible two-dimensional 
boundary layer. We choose the streamwise 
coordinate x along the surface. The coordinate y will 
measure distance from the wall. We define a small 

parameter s = ^Jv / U ref L ref , where v, U ref , and 


L re j are viscosity, reference velocity, and reference 
length, respectively. The streamwise coordinate is 
scaled with Z, ^ • while the vertical coordinate V and 
spanwise coordinate z are scaled with 
/ G ref • Th e following scaling is assumed for 

the velocity disturbances u, v, and w , and the 
pressure p : 


u ~ U ref ’ V ~ sU ref ’ 

w ~ sU re f , p~e 1 pU ref 


( 1 ) 


This scaling of the linearized Navier-Stokes 
equations and neglecting of the curvature effects lead 
to the governing equations for Gortler instability with 
the Gortler number equal to zero. We look for a 
periodic solution in the spanwise direction with the 
corresponding wave number p . The governing 
equations for the amplitude functions can be written 
in dimensionless form as follows: 5 ’ 6 


where u(x,y^) and F(x, y) are the streamwise and 

normal velocity components of the base flow, 
respectively (in addition, the latter is divided by s ). 
The following boundary conditions are applied to the 
solutions: 


y = 0: 

s 

II 

< 

II 

II 

o 

(6a) 

T co : 

u, w, p — > 0 

(6b) 


The equations (2) to (5) can be solved subject to 
boundary conditions (6a) and (6b) with prescribed 
initial velocity perturbations at x = x tl . 

The governing equations can be recast as 
follows: 


(Af) x =B 0 f + B l f y+ B 2 f yy ( 7 ) 

where A, B 0 , B^andB-, are 4x4 matrices (one 
can find them in the Appendix; see also Ref. 5) and 

f = (u,v, w,pY . The superscript “ T ” stands for 

“transposed,” and the subscripts “x” and “y” denote 
differentiation with respect to x and y, respectively’. 

Optimization of Energy Growth 

The authors of Refs. 5 and 6 employed an iterative 
procedure to find the optimal disturbances in terms of 
the maximum of the energy growth ratio 

G = E out /E in’ where E in andE out stand for the 

input and output energy norms. Andersson et al. 5 
used the same definitions of E jn and E out as for the 
disturbance energy 


du dv 

1 -+ Pw = 0 

dx dy 


( 2 ) 


d / x du dU d u 2 

— {Uu) + V ! -r— + v— = — -p 2u (3) 

OX dy dy Sy~ 

— (uV + vU)* — (2Fv) + 

dx dy 

2 ( 4 ) 
„ dp d v i 

PVw-\ = P v 

dy dy 2 

3 / x 8 , \ d 2 w 9 

(Uwj- 1 ( Vwj~Pp = — —~p~w ( 5 ) 


dx 


dy 


dy- 


E = I 
0 

whereas Luchini 6 employed the knowledge that the 
optimal disturbances are represented by streamwise 
vortices with corresponding output as the streamwise 
velocity streaks, 


2 2 2 2 2 
U + £ V + S W 


dy 


( 8 ) 


r\ y max / ^ a \ 

E in =e 2 j ( v + w~ ) dy 

(9a) 

o' ’ 


y max 2 


E out = 1 u dy 

(9b) 


0 
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y max 2 

1 u dy 

G _ ^out_ _ s ~2 _ 0 (9c) 

E- y max / ~ \ 


As was shown in Ref. 5, the two definitions of 
the optimal disturbances lead to the same results at 
Reynolds numbers of 10 4 and higher. Because the 
iteration procedure based on the optimization of ratio 
(9c) provides significant simplification, we adopt it 
for the following analysis. Because Eqs. (2) to (5) are 

2 

independent of s , the value of s G is invariant with 
respect to the Reynolds number. 

Numerical results 

Falkner-Skan Base Flow 

We consider a Falkner-Skan family of boundary- 
layer profiles with free-stream velocity distribution 

U ’ = Cx m and corresponding Hartree parameter 

P H = 2m / (m + 1) . For convenience, we have used 


0.2 and 1.0, respectively. The Reynolds number Re^ 
in Fig. 1 and what follows is defined as U eL L / v . 
One can see that an unfavorable pressure gradient 
{Ph < 0 j leads to an increase in the energy growth 

while a favorable pressure gradient ( P H > 0] leads 

to suppression of the transient growth mechanism. 
The latter is consistent with results obtained within 
the scope of parallel flow approximation. 5 

Figures 2 and 3 show similar results, but the 

starting points are x in / L =0.4 and 0.6, respectively. 
A comparison of Figs. 1 to 3 indicates that there is a 
spanwise wave number, pH re j- = 0.5 -0.6, and a 

starting point, x jn / L , that maximize the energy 

growth. These parameters correspond to optimal 
spanwise spacing and streamwise placing of 
perturbators for maximum flow control effectiveness. 
Figure 4 demonstrates the energy ratio versus the 

downstream coordinate x/L at PH re f = 0.5 and 
x in IL =0.2. 


the velocity scale U re j- = U eL = CP" and the length 
scale L re j ■ = LI (m + 1) . The latter allowed the use of 
the conventional scaling of boundary-layer solutions 
with H ref = ^vLl(m + \)U eL = ]vL ref !U ref . 

Figure 1 shows the scaled energy ratio versus 
spanwise wave number P for three Flartree 

parameters, P H = -0.1, 0.0, and 0.1. The starting and 

the ending points, x in / L and x out / L , are equal to 


Example of LPT Conditions 

Volino 13 simulated low-pressure turbine (LPT) 
airfoil conditions in a low-speed wind tunnel. The 
test section was designed as a passage between two 
airfoils. The local free-stream velocity at the 
favorable pressure-gradient region was closely 
approximated by the following equation: 


U„ 


= 1.48 


U 


exit 


( ^ 0.214 


\ L sJ 


( 10 ) 



Figure 1. Effects of the spanwise wave number p 



and the Hartree parameter p H on transient growth Figure 2 , The same as Fig. 1, with starting point 
(starting point x in / L =0.2 ). x jn / L = 0.4. 


NASA/TM— 2003-212288 


4 

American Institute of Aeronautics and Astronautics 






P H L, 


Figure 3. The same as Fig. 1, with starting point 
Xi n IL = 0 . 6 . 



Figure 4. Optimal energy ratio versus the 
downstream coordinate at three Flartree parameters 
(x ln /L = 0 . 2 , (m L = 0 . 5 ). 

where L s is the suction surface length and U exit is 

the nominal exit free-stream velocity based on the 
inviscid solution. The distribution (10) corresponds to 
the Falkner-Skan flow with the Flartree parameter 
Ph= 0.353. 

Figure 5 demonstrates the energy ratio scaled 
with the Reynolds number R e ex! - f = U exit L s t v versus 
the spanwise wave number scaled with 

H Ls = yjvL s / U exit . The ending point was 

prescribed at x out IL s = 0.444 while the starting points 
varied from 0.111 to 0.289. The streamwise velocity 
perturbation at x out /L s = 0.444, Xi n /L s = 0.111, and 
PHls = 0.925 is shown in Fig. 6, and the 
corresponding optimal profiles of v and w are 
presented in Fig. 7. 

The results indicate that we are dealing with a 
very strong favorable pressure gradient that 


Figure 5. Effects of the spanwise wave number p and 
the starting point Xi n /L on transient growth at 
conditions of the experiment in Ref. 13 (Ph = 0.353, 
XoutlL s = 0.444). 



Figure 6 The streamwise velocity perturbation at the 
ending point x out /L s = 0.444. The parameters 
correspond to the experimental conditions in Ref. 13 
( x in!L s = 0.111, PHls = 0.925). 

suppresses the transient growth mechanism. For an 
example of a typical LPT cruise Reynolds number of 
50,000, the transient growth will provide an energy 
amplification of less than 50. This is a relatively 
small number. If we take into account that in practice 
the perturber will not produce the optimal inflow 
field, the real amplification will be even of a smaller 
value. For example, in Blasius boundary layer, the 
theory predicts amplification of 250 at the same 
Reynolds number of 50,000 (Ref. 9). Correlation 
between the transient growth factor and transition has 
not been established yet; therefore the effectiveness 
of the transient growth mechanism in preventing flow 
separation cannot be assessed quantitatively 
currently. 


NASA/TM— 2003-212288 


5 

American Institute of Aeronautics and Astronautics 








Figure 7. The optimal velocity perturbations at 
XiJLs = 0.111 and PHl s = 0.925 corresponding to the 
streamwise velocity perturbation at x 0Ut /L s = 0.444 
shown in Fig. 6. 



: , 1 1 1 1 " ~ y >-i 

0.0 0.5 1.0 1.5 2.0 2.5 

x/L s 


Figure 8. Effect of the temperature factor on energy 
growth at the experimental conditions of Ref. 13. 


There is a possibility of enhancing the transient 
growth mechanism by means of wall cooling. The 
effect of wall cooling was investigated by Tumin and 
Reshotko 9 within the scope of a parallel flow 
approximation. In order to estimate possible increases 
of the energy ratio on a cold wall at a high favorable 
pressure gradient, we utilize the method of Ref. 9 for 
a compressible flow with local Mach number 0.5 and 
Hartree parameter 0.353. The results are shown in 
Fig. 8. One can see that cooling of the wall might 
provide a tenfold increase in the energy ratio. 

Summary 

The results for the transient growth phenomenon 
within the scope of the linearized boundary-layer 
equations in the presence of a streamwise pressure 
gradient are consistent with previous results obtained 


within the scope of the parallel flow approximation 
and linearized Navier-Stokes equations. 9 A favorable 
pressure gradient decreases the non-modal growth 
while an unfavorable pressure gradient leads to an 
increase of the amplification. 

The example of a Falkner-Skan flow with a 

Flartree parameter = 0.353 corresponds to 

experimental data 13 and simulates the flow over a 
low-pressure turbine airfoil upstream of the 
separation point. At this pressure gradient, the 
transient growth mechanism is suppressed, and the 
energy amplification at low Reynolds number has a 
small value. The theory of the transient growth 
mechanism predicts that it is possible to enhance the 
energy growth by means of wall cooling. The 
example within the scope of the parallel flow theory 9 
demonstrates that cooling of the wall might provide a 
tenfold increase in the energy ratio. Future 
experiments on boundary layer tripping accompanied 
by wall cooling will contribute to our understanding 
of the bypass transition mechanism. 

The method predicts that there is an optimal 
spacing between perturbers and their optimal location 
from the leading edge. The latter results can be 
utilized in future experiments with tripping of the 
boundary layer over the LPT airfoil. These type of 
experiment are planned to be carried out in facilities 
at NASA Glenn Research Center, 14 at the US Naval 
Academy, 13 and at the University of Notre Dame. 13 

Consideration of the optimal velocity 
perturbations in Fig. 7 indicates that they are 
distributed across the boundary layer. This means 
that an array of generators localized on the wall will 
not provide excitation of the optimal disturbances. 
Therefore, the question of realizability of the optimal 
disturbances arises. For example, one can solve the 
receptivity problem for an array of generators on the 
wall and find generator shapes (or other parameters), 
that provide velocity disturbance profiles closest to 
the optimal ones. Another option is to solve the 
receptivity problem for distributed generators 
upstream of the starting point, Xj n , and to find a 
generators distribution, that leads to the optimal 
disturbances. The next option is to design a 
disturbance generator that directly affects the flow 
inside the boundary layer instead of perturbing the 
near-wall region only. For example, it might be a 
focused laser beam projected from the wall, to where 
it could be delivered by a fiber-optic system. 16 These 
fundamental issues should be addressed in future 
research programs on the application of bypass 
transition mechanisms to separation flow control at 
low Reynolds numbers. 
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Appendix 

For purpose of consistency of the formulation, we 
repeat the main features of the numerical scheme 
described in Ref. 5. Matrices A, B 0 , Bj,andB 9 in 
Eq. (7) are as follows: 



( 0 
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0^ 



- p 2 

~ U y 
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0 - 

0 - 

'2 V y -p- 
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-pv 
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l 0 
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~ V y 
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p) 


6 0 

-1 
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O 



6 o 
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0^ 

-V 
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0 -1 
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0 

1 




f 1 0 
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The scheme utilized for Eq. (7) in the streamwise 
direction is 


(Af) 1 -(Af)° = 


Ar 


( b 0 /) 1 +( b 1 /^) 1 + ( b 2 4 ; ) 1 


(Afr- 2 (Af)- + i(Afr= 


2 

Ar 


2 
n + 1 


(Bo/)" + I +(B Jy) n+ +(B 2 f yy ) 




, n > 1 


where n stands for the step number along the 
coordinate x. At each streamwise position, the one- 
dimensional boundary-value problem is solved using 
a spectral collocation method based on Chebyshev 
polynomials. Usually, we used 100 intervals along 
the coordinate x and 100 Chebyshev polynomials for 
the solution approximation. 

To find the optimal perturbations corresponding 
to the maximum of the energy ratio in (9c), the 
forward solution of Eq. (7) is accompanied by the 
backward solution of the adjoint problem. 5 ' 6 In the 
present work, the adjoint system was discretized, i.e., 
the discretized adjoint equations were employed. 
Another approach is to utilize the adjoint form of the 
discretized forward equations (see discussion in 
Ref. 5). An arrangement of the iterations is described 
elsewhere. 5 ' 6 Usually, 2 to 3 iterations were enough to 
achieve convergence at the 0.1 percent level. 
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